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Abstract 

We discuss the implications of multi-brane constructions involving combinations of pos- 
itive and negative tension brane and show how anomalously light KK states emerge when 
negative tension "— " branes are sandwiched between branes. We present a detailed 

study of a " H V" brane assignment which interpolates between two models that have 

been previously proposed in which gravity is modified at large scales due to the anomalously 
light states. We show that it has the peculiar characteristic that gravity changes from four 
dimensional (4D) to 5D at large distances and returns to 4D at even larger scales. We also 
consider a crystalline universe which leads to a similar structure for gravity. The problems 
associated with intermediate negative tension branes are discussed and a possible resolution 
suggested. 
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1 Introduction 



The conjecture that we may live on a brane in a space-time with more than four dimen- 
sions is quite old but has been the subject of renewed interest in recent years with the 
realisation that such structures are common in string thories. The models of Antoniadis, 
Arkani-Hamed, Dimopoulos, Dvali (AADD) 0] and of Randall, Sundrum (RS) || [|] have 
the common characteristic that the Standard Model (SM) fields are localized on a 3-brane 
and gravity propagates in all the space-time dimensions (the "bulk"). They can provide 
novel explanations to the Planck hierarchy problem based either on the size of the bulk 
volume (AADD), or on an exponential warp factor that enters the metric (RS). The RS 
construction in particular is very attractive because all the parameters of the model may 
be of the same magnitude while still generating a very large hierarchy. 

The RS model consists of two 3-branes of opposite tension sitting at the fixed points of 
an S 1 /Z 2 orbifold with AdS$ bulk geometry. Gravity is localized on the brane of positive 
tension ("+" brane) whereas an exponential "warp" factor generates a scale hierarchy on 
the brane of negative tension ("— " brane). There are two possible physical interpretations 
of this model. Placing the Standard Model (SM) fields on the "—" brane provides an 
explanation of the Planck hierarchy problem. Placing the SM fields on the "+" brane [Q 
allows the "— " brane to be put at infinity, i.e. decompactify the orbifold. In this case, 
however, gravity is localised on the brane so the world still appears four dimensional 
(see also |J). As far as the the first possibility is concerned, there has been a lot of 
discussion || whether we can have acceptable theory of gravity and cosmology if we live 
on a "— " brane. To avoid this problem Lykken and Randall (LR) proposed a " + H — " 
brane configuration in which the SM fields reside on the intermediate positive tension brane 
(which has vanishingly small tension) and which still has a warp factor explanation of the 
mass hierarchy. 

In reference p) we considered an alternative " H [-"configuration, again with the SM 

fields on a positive tension brane and a warp factor. This model (KMPRS) has a novel 
phenomenology || quite distinct from the RS and LR models. This is mainly due to the 
intermediate "— " brane leading to an anomously light first Kaluza-Klein (KK) as discussed 
below. This feature gives rise to the exotic possibility of "bigravity" , in which the first KK 
state has a Compton wavelength of the order of 10 26 cm corresponding to 1% of the size of 
the observable universe, while, at the same time, all the remaining states of the KK tower 
have Compton wavelengths below 1mm. In such a scenario the gravitational attraction as 
we feel it is the net effect of the exchange of the ordinary 4D graviton and the ultralight 
KK state. This leads to the prediction of modifications of gravity not only at the millimeter 
scale (due to the higher KK modes) but also the ultralarge scale (of O(10 26 cm)!) Indeed, 
gravity becomes weaker at ultralarge distances due to the Yukawa suppression of the first 
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KK mode. 

Independently and again in the context of the RS scenario, Gregory, Rubakov, Sibiryakov 
(GRS) have recently suggested [0 a construction in which gravity is also modified at both 



small and at ultralarge scales. The GRS model modifies the decompactified RS model Q by 
adding a "— " brane of half the tension of the "+" brane and requiring flat space to the right 
of the new "—" brane. This model does not have a normalizable 4D graviton but generates 
4D gravity at intermediate distances due to a resonance- like behaviour of the wavefunctions 
of the KK states continuum. Gravity in this picture appears to be "quasi-localized" and 
for ultra large scales becomes five dimensional. In this scenario gravity is also modified at 
short scales (but this time much shorter than the millimeter scale). 

Although the KMPRS and GRS models look quite different, they share the key ele- 
ment of freely moving "— " branes which generate the tunneling effects responsible for the 



anomalously light states. It was shown in [|Tl|] that these two models are the limiting cases 

of a more general " H h" multi-brane model that interpolates between the "bigravity" 

KMPRS model and the "quasi-localized" gravity GRS model. This model can be readily 
derived from the KMPRS model by cutting in half the intermediate "— " brane and consid- 



ering flat space between the two new "— " branes. Moreover, in [] 1 it was conjectured that 
every construction having "— " branes between "+" branes is capable of giving a "multi- 
gravity" scenario where gravity is generated not only by the 4D graviton but also by a 
number (finite or infinite) of KK states. The crystal universe considered by Jl2j (see also 
|i~3| ) is an interesting example^. Let us mention here that recently Gorsky and Selivanov 
Hi~4| discussed models where the effect of negative tension branes can be mimicked by a 
constant four form field and the presence of the three brane junctions. It will be interesting 
to see if in these models one has multigravity. 

The purpose of this paper is threefold. We first we will consider the " + H — " LR 
construction but with the intermediate brane carrying non-negligible tension. We show 
that this model does not have an anomously light KK state and hence cannot generate 
a "multigravity " scenario. Analysis of the reason for this suggests that branes not 
separated by "— " branes do not lead to anomalously light modes. Secondly, we examine in 

detail the " H \-" "multigravity" model and show exactly how it interpolates between 

the KMPRS and the GRS model. Gravity in this case has the quite peculiar characteristic 
that it changes from 4D to 5D at ultra-large scales and then again to 4D at even larger 
scales. We also analyse the infinite crystalline universe and show how the band structure 
of the KK modes gives rise to "multigravity" . Finally we discuss the role of the radion and 
the question whether these theories give rise to viable theories. 



3 We disagree with the statements made in [Q about the absence of a zero mode band. 
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2 The three-brane " + +-" Model 

The " + H — " model (see Fig.[l]) consists of three parallel 3-branes in AdS§ spacetime with 
orbifold topology, two of which are located at the orbifold fixed points Lq = and L 2 
while the third one is moving at distance L\ in between. In order to get 4D flat space 
with this configuration, it turns out that the AdS 5 space must have different cosmological 
constants Ai and A2 between the first - second and the second - third brane respectively 
with I A.2 1 > I Ai| (see [|nj for constructions of different bulk cosmological constants). The 
action of the above setup (if we neglect the matter contribution on the branes in order to 
find a suitable vacuum solution) is given by: 

[ d i xV t \[^) (1) 



where G^l is the induced metric on the branes, are their tensions and M the 5D fun- 
damental scale. We consider as in || the vacuum metric ansatz that respects 4D Poincare 
invariance: 

ds 2 = e- 2a{y) r J „ u dx^dx v + dy 2 (2) 




Figure 1: The brane locations in the three-brane " + H — " model. The bulk curvature 
between the "+" branes is k\ and between the "+" and "— " brane is k 2 . 

The above ansatz substituted in the Einstein equations requires that a(y) satisfies the 
differential equations: 

*" = Yl^h-Ay-Li) ^d (a') 2 =(g ' Z/G !° r ' Ll r ] 1 (3) 
^ 12M 3 vy ; 1*2 »ye Li,La v; 

i K 



s 

with A(y) 



L 2 



d 4 x I dyV^G[-A(y) + 2M 3 R] 

J-L 2 

Ai ,ye[0,L!] 
A 2 ,ye[L u L 2 } 
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where k\ = y 2A ^ 3 and k 2 = y 2 ^z are effectively the bulk curvatures in the two regions 
between the branes. The solution of these equations for y > it is given by: 

a[y) -\ k 2 (y - + hL, ,ye[L u L2\ [V 



Einstein's equations impose the following conditions on the brane tensions : 

k — k 

V = 2AM 3 k u Vx = 2AM 3 2 1 , V 2 = -2AM 3 k 2 (5) 



If we consider massless fluctuations of this vacuum metric and then integrate over the 
5-th dimension, we find the 4D Planck mass is given by: 



M|j = M 3 



— (l - g-Sfci^A + J_ e 2(fc 2 -fci)£i f e -2k 2 Ll _ e -2k 2 L 2 \ 

k\ k 2 ^ 



(6) 



The above formula tells us that for large enough kh\ and kL 2 the four mass scales Mpi, 
M, ki and k 2 can be taken to be of the same order. Thus we take hi, k 2 ~ O(M) in order 
not to introduce a new hierarchy, with the additional restriction k\ < k 2 < M so that the 
bulk curvature is small compared to the 5D Planck scale and we can trust the solution. 
Furthermore, if we put matter on the second brane all the physical masses m on it will be 
related to the mass parameters mo of the fundamental 5D theory by the conformal "warp" 
factor 

m = e~ CT(Ll) m = e~ klLl m (7) 

This allows us to put the SM states on the intermediate "+" brane, solving the Planck 
hierarchy problem by choosing e~ kLl to be of O(10~ 15 ), i.e L\ w 35fc _1 . 

The KK spectrum can be as usual found by considering the linear "massive" fluctuations 
of the metric. We use the following convention^]: 



ds 2 



e ~ Mv) V,u + j^K u (x,y) 



dx^dx u + dy 1 



We expand the field h IMV (x, y) in terms of the 4D graviton zero mode and the KK states 
plane waves h$(x): 

00 

h, u (x,y) = Y,h ( ;J(x)¥ n) (y) (9) 

n=0 

where (d K d K — m^) h$ = after fixing the gauge as cPh^ = h & = 0. The wavefunction 
ty( n \y) obeys a second order differential equation and carries all the information about the 



6 Hcre we have ignored the radion modes that could be used to stabilize the brane positions L\ and L^. 



For discussion and possible stabilization mechanisms see |16l 
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effect of the non-factorizable geometry on the graviton and the KK states. After a change 
of coordinates and a redefinition of the wavefunction the problem reduces to the solution 
of an ordinary Schrodinger equation: 



{-^ + 7(z)}*wW = !|*W| 



(10) 



where the potential V(z) for z > has the form 

15 



V(z) 



3 

~2^y 



[k 2 (9(z) - 9(z - z x )) + k 2 2 (6(z - Zl ) - 6(z - z 2 ))] 



ki6(z) + ^ 2 - h) 5(z - zi) - k 2 5(z - z 2 ) 



The new coordinates and wavefunction in the above equations are defined by: 

,ye[o,Li] 



z = 



kl 

£ k 2 (y-L 1 ) + k 1 L 1 e k\ L l —\ e k l L l 

r TT. 



k2 ' kl k 2 ,y£[Li,L 2 ] 

¥ n) {z) = ¥ n) {y)e a/2 
with the symmetric z for y < and the function g(z) as: 



9W 



kiz +1 , z G [0, zi] 

k 2 (z - zi) + k x z x + 1 ,ze[zi,z 2 ] 



(12) 
(13) 

(14) 



where z x = z(Li) and z 2 = z(L 2 ). The change of coordinates has been chosen so that there 
are no first derivative terms in the differential equation. Furthermore, in this coordinate 
system it can be easily seen that the vacuum metric takes the conformaly flat form: 



ds 2 



{ri^ydx IJ 'dx lJ + dz 2 ) 



(15) 



The potential ( |Tl~l) always gives rise to a massless graviton zero mode which reflects the 
fact that Lorentz invariance is preserved in 4D spacetime. Its wavefunction is given by: 



A 



with normalization factor A determined by the requirement 



dz 



2 12 



#(°)( 



(16) 



1, chosen 



so that we get the standard form of the Fierz-Pauli Lagrangian (the same holds for the 
normalization of all the other states). 
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For the massive KK modes the solution is given in terms of Bessel functions. For y lying 
in the regions A = [0, L±] and B = [Li, L 2 ], we have: 



A 
B 



N„ 



£(£) 
ki 

fc 2 



(17) 



The boundary conditions (one for the continuity of the wavefunction at Z\ and three 
for the discontinuity of its first derivative at 0, zi, z 2 ) result in a 4 x 4 homogeneous linear 
system which, in order to have a non-trivial solution, should have a vanishing determinant: 







Y i (f ) J * (f ) 

Y, (f 2 g{z 2 ) 

^(f^l)) J 2 (f^l)) -v^(i^l)) - 

(The subscript n on the masses m n has been suppressed.) 







(18) 



2.1 Masses and Couplings 

The above quantization condition determines the mass spectrum of the model. The pa- 
rameters we have are k±, k 2 , L±, L 2 with the restriction k\ < k 2 < M and ki ~ k 2 so 
that we don't introduce a new hierarchy. It is more convenient to introduce the parameters 
x = k 2 (L 2 — Li), w = e~ klLl and work instead with the set ki, k 2 , x, w. From now on 
we will assume that w < 1 (to ~ (9(10 -15 ) as is required if the model is to provide an 
explanation of the hierarchy problem). 

For the region x > 1 it is straightforward to show analytically that all the masses of the 
KK tower scale in the same way as x is varied: 

m n = i n k 2 we~ x (19) 

where £ n is the n-th root of J\{x). This should be compared with the KMPRS "H h"model 

in which m\ oc kwe~ 2x and m n+ i oc kwe~ x , where x is the separation between the "— " and 
the second "+" brane. This significant difference can be explained by the fact that in the 
"-H — "case the negative tension brane creates a potential barrier between the two attractive 
potentials created by the positive tension branes. As a result the wave function in the region 
of the "— " brane is small due to the tunneling effect. The two attractive potentials support 
two bound states, one the graviton and the other the first KK mode. The mass difference 
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between the two is determined by the wave function in the neighbourhood of the "— " brane 
and is thus very small. On the other hand the wave function between the two branes 
in the " + H — " configuration is not suppressed by the need to tunnel and hence the mass 
difference between the zero mode and the first KK mode is also not suppressed. This has 
as result the two "+" branes behave approximately as one. This becomes even more clear 
when x ^> f where the model resembles the " H — " RS construction. Indeed, in this limit 
the mass spectrum becomes m n = £, n k 2 e~ k2L2 which is exactly that of the " H — " RS model 
with orbifold size L 2 and bulk curvature k 2 . 

In the region x < 1 the relation of eq(^) breaks down. As reduce x the second "+" 
comes closer and closer to the "—" brane and in the limit x — > {i.e. L 2 = L\) the combined 
brane system behaves as a single "—" brane, reducing to the " H — " RS model. In this limit 
the spectrum is given by: 

m n = inkiw (20) 

which is just the one of the " H — " RS model. In the region ^ x < 1 the mass spectrum 
interpolates between the relations (|19"D and ((20|). 

The fact that there is nothing special about the first KK mode is true also for its 
coupling on the second "+" brane. The interaction of the KK states on the second "+" 
brane is given by: 



C int = ^2a n h^J(x)T f , L ,(x) , with a n = 



n>0 



9(zi) 
M 



3/2 



\]>( n )( 



(21) 



In the RS limit (x = 0) all the states of the KK tower have equal coupling given by: 

1 



22 



pi 



As we increase x, the lower a state is in the tower, the more strongly it couples, i.e. 
a i > a 2 > a 3 > ' ' ' (with di < (uiMpi) -1 ) and tends to a constant value for high enough 
levels. At some point this behaviour changes, the levels cross and for x > 1 the situation 
is reversed and the lower a state is in the tower, the more weakly it couples, i.e. a\ < 
a 2 < as < ■ ■ ■. In this region it is possible to obtain a simple analytical expression for the 
couplings: 

n "j 2 y UJ wM Pl e {Z6) 

Here, we also see that the first KK state scales in exactly the same way as the remaining 
states in the tower with respect to x, a behaviour quite different to that in the KMPRS 



- 3:c 



i.e. 



model in which the coupling is x-independent. Furthermore, the coupling falls as e 
much faster than e~ x as one would naively expect. This can be explained by looking at the 
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origin of the x-dependence of the wavefunction. For increasing x the normalization volume 
coming from the region between L\ and L2 dominates and the normalization constant in 
(fl7|) scales as N n oc e~ 3x . This rapid decrease is not compensated by the increase of the 
value of the remaining wavefunction (which from is approximately constant). Thus, 
although the two "+" branes in the large x limit behave as one as far as the mass spectrum 
is concerned, their separation actually makes the coupling of the KK modes very different. 



3 The four-brane " + +" Model 

In this section we discuss the four-brane model in order to clarify the relation between the 
"bigravity" KMPRS model || with the GRS model |nj. In particular we wish to explore 



and compare the modification of gravity at large scales predicted by each model. 

In the case of the KMPRS " H V" model, in the limit of very large x, gravity re- 
sults from the net effect of both the massless graviton and the ultralight first KK state. 
The modifications of gravity at very large distances come from the fact that the Yukawa 
type suppression of the gravitational potential coming from the KK state turns on at the 
Compton wavelength of the state. On the other hand, the GRS model has a continuous 
spectrum with no normalizable zero mode. However, the values of the KK states wavefunc- 



tions on the "+" brane have a "resonance" -like behaviour |17] which give rise to 4D gravity 
at distances smaller than the Compton wavelength of its width. Beyond this scale gravity 
becomes five- dimensional. 

The four-brane GRS configuration can be obtained from the KMPRS model by "cutting" 
the "— " brane in half, i.e. instead of having one "— " brane one can take two "—" branes of 
half the tension of the original one (" — 1/2" branes), having flat spacetime between them 
(see Fig.(|2|). Finally if the second "+" brane is taken to infinity together with one of the 
" — 1/2" branes we shall get precisely the GRS picture. 

Let us discuss the four-brane "H h" model in more detail. It consists of 5D spacetime 

with orbifold topology with four parallel 3-branes located at L = 0, L 1; L 2 and L 3 , where 
L and L 3 are the orbifold fixed points (see Fig.(^)). The bulk cosmo logical constant A 
is negative (i.e. AdS 5 spacetime) between the branes with opposite tension and zero (i.e. 
flat spacetime) between the two " — 1/2" branes. The model has four parameters namely 
Li, L2 and L3 and A. For our present purposes we consider the symmetric configuration, 
leaving 3 parameters, I, Z_ and A where I = L\ = L 3 — L 2 and Z_ = L 2 — L±. In the absence 
of matter the model is described by eq(|l]) with 

A(y) = I ^ [^1 ? ^2] m 
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-1/2 -1/2 

Figure 2: " -\ h" configuration with scale equivalent "+" branes. The distance between 

the "+" and " — 1/2" branes is I = L\ = L 3 — L 2 while the distance between the " — 1/2" 
branes is £_ = L 2 — L\. The curvature of the bulk between the "+" and " — 1/2" branes is 
k. 



By considering the ansatz eq@ the "warp" function cr(y) must satisfy: 



,yE [Li, L 2 ] 

k 2 ,ye[0,Li]\J[L 2 ,L 3 ] 



(25) 



where k = J 24j ^ 3 is a measure of the bulk curvature and we take V = V 3 = —2V\ 
— 2V 2 = V. The solution for y > is: 



ky ,ye [0, L x ] 

a{y) = { kLi ,y e [L 1 ,L 2 \ 

kL x + k(L 2 -y) ,ye[L 2 ,L 3 ] 



(26) 



Furthermore, 4D Poincare invariance requires the fine tuned relation: 

V = 



A 

T 



(27) 



In order to determine the mass spectrum and the couplings of the KK modes we consider 
linear "massive" metric fluctuations as in eq([|). Following the same procedure we find that 
the function iff( n '(z) obeys a Schrodinger-like equation with potential V(z) of the form: 



V(z) 



15k 2 



3k 



[9(z) - 6{z - Zl ) + 9(z - z 2 ) - 9(z - z 3 )} 
K z ) - 2 S ( Z ~ Zl ) ~~ \^ z ~ Z2 ) + S ( z ~~ 



(28) 
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The conformal coordinates now are given by: 



z = 



e k y-i 

k 

(y - l)e kl + ^ 

j_g2fc/+/c/_ g— ky _|_ ^ gfci _|_ 2^kl 1 



,ye[o,Li] 

,y e [L X ,L 2 ] 
,y e [L 2 ,L 3 ] 



(29) 



hz + 1 , z e [0, zi] 

.</(:) = <; fc«i + i ,ze[z u z2] (30) 

fe(z 2 - z) + kzx + 1 , -2 € [-22, -23] 

where zi = -2(^1), -22 = -2(^2) and Z3 = z(L 3 ). 

The potential f^8|) again gives rise to a massless graviton zero mode whose wavefunction 
is given by fll6|) with the same normalization convention. Note, however, that in the limit 
/_ — > 00 this mode becomes non-normalizable (GRS case). The solution of the Schrodinger 
equation for the massive KK modes is: 




>/¥t Y 2 (^g(z))+AJ 2 (^g(z))} ' 
Bi cos(m n z) + i?2 sin(m n z) / 
^ [CiVa (^)) + C 2 J a (^g(z))] 



(31) 



where A = [0, z%], B = [2:1,-22], and C = [-22,-23]- We observe that the solution in the first 

and third interval has the same form as in the " H V" model. The new feature is the 

second region (flat spacetime). The coefficients that appear in the solution are determined 
by imposing the boundary conditions and normalizing the wavefunction. 

The boundary conditions (two for the continuity of the wavefunction at z x , z 2 and four for 
the discontinuity of its first derivative at 0, z\, z 2 and z 3 ) result in a 6 x 6 homogeneous linear 
system which, in order to have a non-trivial solution, should have vanishing determinant. 
It is readily reduced to a 4 x 4 set of equations leading to the quantization condition: 



n(f) 






- cos(mzi) 
sin(mzi) 

- sin(?7iZ2) 
cos(mz2) 



- siii(mzi) 
cos(mzi) 
cos(mz2) 
sin(mz2) 



* Oaf ) 

-Y 2 (g 2 f) 






J l(f3x) 

Mo*?) 



J 2 (.g 2 f) 



(32) 



with gi = g(zi), g 2 = g{z 2 ) and g$ = g{z%). Here we have suppressed the subscript n on 
the masses m n . 
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3.1 The Mass Spectrum 



The above quantization condition provides the mass spectrum of the model. It is convenient 
to introduce two dimensionless parameters, x = kl and x_ = kl- (c.f. Fig.([|)) and we work 
from now on with the set of parameters x, X- and k. The mass spectrum depends crucially 
on the distance x_. We must recover the KMPRS spectrum in the limit x_ — > 0, and the 
GRS spectrum in the limit x_ — > oo. From the quantization condition, eq(|32"D it is easy to 

verify these features and show how the " H \-" spectrum smoothly interpolates between 

the KMPRS model and the GRS one. It turns out that the structure of the spectrum has 
simple X- and x dependence in three separate regions of the parameter space: 

3.1.1 The three-brane " H V" Region 

For x_ <, 1 we find that the mass spectrum is effectively ^.-independent given by the 
approximate form: 

To! = 2V2ke~ 2x (33) 
m n+1 = Uke~ x n = 1, 2, 3, . . . (34) 

where £21+1 is the {i + l)-th root of J\(x) {i = 0, 1, 2, . . .) and £ 2 i is the i-th root of J2{x) 
(i = 1,2,3, . . .). As expected the mass spectrum is identical to the one of the KMPRS 
model for the trivial warp factor w — 1. The first mass is manifestly singled out from the 
rest of the KK tower and for large x leads to the possibility of bigravity. 

3.1.2 The Saturation Region 

For 1 « i_ <C e 2x we find a simple dependence on x_ given by the approximate analytic 
form: 

-2x 

m 1 = 2k (35) 

to„ +1 = rnik n = l,2,3,... (36) 

X- 

As X- increases the first mass decreases less rapidly than the other levels. 

3.1.3 The GRS Region 

For x_ ^> e 2x the first mass is no longer special and scales with respect to both x and x_ 
in the same way as the remaining tower: 

m n = nirk n= 1,2,3,... (37) 

X- 
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The mass splittings Am tend to zero as X- — > oo and we obtain the GRS continuum of 

states 

The behaviour of the spectrum is illustrated in Figure 3. 




Figure 3: The behaviour of the mass of the first five KK states in the three regions of simple 
x, X- dependence. The first dot at zero stands for the graviton. 



3.2 Multigravity 

Armed with the details how the spectrum smoothly changes between the KPMRS model 
(x- = 0) and the GRS model (x_ — > oo), we can now discuss the possibilities for modifying 
gravity at large distances. The couplings of the KK states with matter on the left "+" 
brane are readily calculated by the interaction Lagrangian ( pT|) with: 

3/2 

^ (n) (0) (38) 

3.2.1 Bigravity Region 

In the KPMRS limit, x_ - 
the 4D graviton: 



9_W 
M 



0, the first KK mode has constant coupling equal to that of 

1 2M 3 
= Jf (= «o) where Ml = — (39) 
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while the couplings of the rest of the KK tower are exponentially suppressed: 

a n+ i = 4r ~r 6 — n = l,2,3,... (40) 

The gravitational potential is computed by the tree level exchange diagrams of the 4D 
graviton and KK states which in the Newtonian limit is: 

"A p -m n r 

V(') = -XX-T- (41) 

ra=0 

where a n is the coupling fl38|) and n = accounts for the massless graviton. The summation 
stops at some very high level N\ with mass of the order of the cutoff scale ~ M. 

In the "bigravity" scenario, at distances r <C mjf , the first KK state and the 4D graviton 
contribute equally to the gravitational force, i.e. 

V*)«-±(U^)«-°2. (42) 



Ml \r 

where G^ = ■ For distances r > m^ 1 the Yukawa suppression effectively reduces gravity 
to half its strength. Astronomical constraints and the requirement of the observability of 
this effect demand that for k ~ Mpi we should have x in the region 65-70. Moreover, at 
distances r < m^ 1 the Yukawa interactions of the remaining KK states are significant and 
will give rise to a short distance correction. This can be evaluated by using the asymptotic 



expression of the Bessel functions in (40) since we are dealing with large x and summing 
over a very dense spectrum, giving: 

Gm kn m r 



V «W = —2->2Z2k— (43) 

n=2 

At this point we exploit the fact that the spectrum is nearly continuum above m 2 and 
turn the sum to an integral with the first factor in (|50]) being the integration measure, i.e. 

= S ^ m ~ > / dm (this follows from eq(|34]) for the asymptotic values of the Bessel 
roots). Moreover, we can extend the integration to infinity because, due to the exponential 
suppression of the integrand, the integral saturates very quickly and thus the integration 
over the region of very large masses is irrelevant. The resulting potential is now: 

The integration is easily performed and gives: 

We see these short distance corrections are significant only at Planck scale lengths ~ k" 1 . 
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3.2.2 The GRS Region 



In the GRS limit, x_ ^> e 2x , we should reproduce the "resonance" -like behaviour of the 
coupling in the GRS model. In the following we shall see that indeed this is the case and 
we will calculate the first order correction to the GRS potential for the case x_ is large but 
finite. 

For the rest of the section we split the wavefunction (RTJ) in two parts, namely the 
normalization N n and the unnormalized wavefunction ty( n >(z), i.e. $( n >(z) = N n ty( n >(z). 
The former is as usual chosen so that we get a canonically normalized Pauli-Fierz Lagrangian 
for the 4D KK modes h$ and is given by: 



N 



1/2 



f* 1 




2 r z z 




2 


/ dz 


¥ n \z) 


+ dz 


¥ n \z) 




' 




J Z\ 







(46) 



The value of ty( n \z) on the left brane is, for m n <C k: 

16k 3 



*fn)(0) 



(47) 



It is convenient to split the gravitational potential given by the relation ( f4l|) into two 
parts: 



N A 



V(r) 



1 e ^2,^,2 fn\ L e at2,t,2 



M 3 



E 



71=1 



-N 2 ¥f jo) — — 



(4f 



n=N x _ 



As we shall see this separation is useful because the first N x _ states give rise to the long 
distance gravitational potential Vu while the remaining ones will only contribute to the 
short distance corrections V s d- 

• Short Distance Corrections 



We first consider the second term. The normalization constant in this region is computed 



by considering the asymptotic expansions of the Bessel functions with argument 
is easily calculated to be: 



g{z 1 )m n 



. It 



32k 3 g( Zl )x„ 



(49) 



If we combine the above normalization with unnormalized wavefunction (f4"7|), we find 



VUr) 



N A 

E 



kix m n e 



M 3 ^ x-e x 2k 

n=N x 



1 + 



(50) 
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Since we are taking x_ ^> e 2x , the spectrum tends to continuum, i.e. N n — > N(m), 
^(n)(0) — > ty(m), and the sum turns to an integral where the first factor in (|50|) is the 
integration measure, i.e. = Y2^ m ~ > I dm (c.f. eq(j3?l)). Moreover, as before we 

can again extend the integration to infinity. Finally, we expand the fraction involving x_ 
keeping the first term in the power series to obtain the potential: 

V sd (r)^-4-J dm?—£[l-±) (51) 



M 3 J mo r 2k \ x. r 

where m = ke~ x . The integral is easily calculated and the potential reads: 

where we identified Gn = for reasons to be seen later. The second part of the above 
potential is the first correction coming from the fact that x_ is finite. Obviously this 
correction vanishes when x_ — > oo. Note that the above potential gives corrections to the 
Newton's law only at distances comparable to the Planck length scale. 

• Multigravity: 4D and 5D gravity 

We turn now to the more interesting first summation in eq (|48|) in order to show that 
the coupling indeed has the "resonance" -like behaviour for Am — > responsible for 4D 
Newtonian gravity at intermediate distances and the 5D gravitational law for cosmological 
distances. This summation includes the KK states from the graviton zero mode up to the 
iV^ -th level. The normalization constant in this region is computed by considering the 
series expansion of all the Bessel functions involved. It is easily calculated to be: 



^2 _ ^ 2 K 



Ag( Zl ) 4 x. 



(53) 



4 g(zi) i x. 

where T = 4ke~ 3x . If we combine the above normalization with the unnormalized wave- 
function (f47|), we find that the long distance gravitational potential is: 

N x 

~ Trh Ah* p-m„r -| 

(54) 



x ^ ixk 4k 2 



2 _|_ El _|_ 8fc 2 
n i" 4 ' g(^i) 4 z- 



M 3 ^— ' X-e x Trq(zi) 3 

n=0 ° v ' 

Again, since we are taking x_ ^> e 2x , the above sum will turn to an integral with — > 
J dm. Moreover, we can safely extend the integration to infinity since the integral saturates 
very fast for m < T/4 = r~ x <C ke x . If we also expand the fraction in brackets keeping the 
first term in the power series, we find the potential: 

„ , s 1 r° , 4k 2 e~ mr 1 

V ld( r ) - ~ T7Z dm 



M 3 J ng( Zl ) 3 r m 2 + ? 

1 f 00 32k 4 e~ mr 1 

+ Wo dm x^g{z x y — {m 2 + lt ? ^) 
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The first part is the same as in the GRS model potential, whereas the second one is the 
first correction that comes from the fact that X- is still finite though very large. Note that 
the width of the "resonance" scales like e~ 3x , something that is compatible with the scaling 
law of the masses (m n = nnk—), since we are working at the region where x_ ^> e 2x , i.e. 
m n <C nnke~ 3x . The above integrals can be easily calculated in two interesting limits. 

For k~ x <C r <C r c the potential is given approximately by : 

C e 2x 

V ld (r<^r c )^ ^ (1-—) (56) 

r X- 

where we have identified Gn = -^s to obtain the normal 4D Newtonian potential. Note 
that since x_ ^> e 2x , the term is indeed a small correction. 

In the other limit, r ^> r c , the integrand is only significant for values of m for which 
the m 2 term in the denominator of the "Breit-Wigner" can be dropped and the potential 
becomes: 

V ld {r » r c ) ~ (1 _ 2 J—) (57) 

7rr z X- 

The fact that Newtonian gravity has been tested close to the present horizon size require 
that for k ~ Mp\ we should have x > 45-50. 

Finally we note that if we take the x_ — > 00 we recover the GRS result 

lim ^ + __+(r,a;_) = y Gfl s(r) (58) 



• Back to 4D gravity 

As we have just seen, probing larger distances than r c , the 4D gravitational potential 
changes to a 5D one. This is the most significant characteristic of the GRS model. In the 
case that x_ is large compared to e 2x but still finite, there is another distinct region of 
interest, namely r ^> m\ x . This follows from the fact that, in this limit, the spectrum is 
still discrete. For distances larger than of the order of the corresponding wavelength of the 
first KK mode, the contribution to gravity from the KK tower is suppressed and thus the 
zero mode gives the dominant contribution, leading to the 4D Newtonian potential again. 
In this case the strength of the gravitational interaction is a small fraction of the strength 
of the intermediate 4D gravity. More precisely, the contribution of the massless graviton is 
^— suppressed and thus vanishes when x_ — > 00, something that is expected since in this 
limit there is no nomalizable zero mode. The gravitational potential in this case is: 

1 1 n 

V ' D(r) = -IP r N W°> (0) = ~r ~_ (59) 
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Figure 4: The behaviour of the coupling, a(m), in the limit of X- 3> e 2x . Three regions 
of interest are indicated. The region m > mo gives rise to short distance corrections. 
The mi *C in *C m c region gives rise to 4D gravity at intermediate distances and 5D 
gravity at ultra large distances. For distances r ^> vrC[ , the zero mode gives the dominant 
contribution and thus we return to 4D gravity. 

Obviously this 4D region disappears at the limit a;_ — > oo since the spectrum becomes 
continuum and thus the 5D gravity "window" extents to infinity. We should note finally 
that for the values of x that we consider here this final modification of gravity occurs at 
distances well above the present horizon. 



4 The Crystal Universe Model 

The last model we consider is the crystal universe introduced in [12|] (see also (HI). It 
consists of an infinite array of parallel 3-branes in a 5D AdS space with cosmological 
constant A. For simplicity we assume that the branes are equidistant with distance, I, 
between two successive branes. Needless to say, in this case all the "+" branes have warp 
factor w = 1 with respect to a "+" brane sitting at the origin of the 5-th dimension 
coordinate. The metric ansatz that has 4D Poincare invariance is again given by eq(|2|) 
where the a(y) function is constrained by the the Einstein equations to have the sawtooth 
form: 

+oo 

a(y) = kY,(y- 2J0 [20(y - 2jl) - 6(y - (2j - 1)1) - 6(y - (2j + 1)/)] (60) 

j=-oo 
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Figure 5: The Crystal Universe made up of an infinite array of "+" and "— " branes with 
lattice spacing I and bulk curvature k. 



The tensions of successive branes are required to be opposite and equal to ±A/fc, where 
k = 



-A 
24M 3 



is a measure of the curvature of the bulk and M the 5D fundamental scale. 



We consider the general fluctuations around the previous vacuum ansatz of the form 
and decompose as usual the perturbation tensor h fJ/l ,(x,y) as: 



h^(x,y) = / dmh^{m,x)i$/{m,y) 



(61) 



where ^>(m,y) is a complex function that arises from the linear combination of the inde- 
pendent graviton polarizations. This leads to the Schrodinger equation fllOD for the wave 
function ^(m, z) = ^(m,y)e a / 2 with potential: 



V(z) 



15k 2 



3k 



8[g(z)] 2 2g(z) 



Z - Zi 



(62) 



j=-oo 



The z coordinates have been defined as usual to be the ones that make the background 

kl i 

metric conformally flat. In these coordinates the branes sit at the points Zj = j e ,~ = jzi 
and the function g(z) is: 



g(z) = l + k }^ (z - z 2j ) [26 (z - z 2 j) - 9{z - z 2j -i) - 9{z - z 2j+ i)] 



(63) 



J=-oo 



The solution of the Schrodinger equation for the wavefunctions for two adjacent cells is 
given [ 



121 in terms of Hankel functions: 



y(m,z) = N n 



i 

B 

e 2iqzi 
,2iqzi ft 




(64) 



in the regions z G [0, zj\, z £ [zi,2zi], z £ [2^,3^] and z £ [3z/,4z ; ] respectively Q. In 
7 We use as [O] Hankel instead of real Bessel functions in order to encode the gravitons phase rotation. 
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the above expression N m is an overall normalization constant, q is the Bloch wave quasi- 
momentum and the constants A,B,C are given by: 

A = (hjh 2 + hjhpht - (hp4 ~ h\h- 2 )e 2 ^hj - 2h\hjh 2 
{h\h 2 + K{h\)h 2 + {h^hj - h\h 2 )e 2i ^h 2 -2hyh 2 h^ 
B = hfh 2 + h^h+ + 2h^h 2 A ^ 

h^h 2 — h^h 2 
c = 2hjh+ + (hjh 2 + h- x hj)A 
h\h 2 — h^h 2 

where = (^g(zi)) and = (^). The above coefficients of the Hankel functions 
were determined by the boundary conditions, i.e. continuity of the wavefunction at zi, 2zi, 
3zi, discontinuity of its first derivative at zi, 3zi and the Bloch wave conditions relating 
the wavefunction at the edges of each cell. The last remaining boundary condition of the 
discontinuity of the first derivative of the wavefunction at 2zi gives us the band equation 
connecting q and m. We disagree with the relation given in [|T^] and instead we find: 

cos(2qzi) = : — r— = f(m) (68) 

U22/1 - jm){j2yi -jm) 

where again j n = J n (fg(zi)), y n = Y n (fg(zi)) and j n = J n (f ), y n = Y n (f ). This 
dispersion relation always gives a band at zero as is to be expected intuitively. Defining 
the parameter x = kl, we see that for x ^ 10 we may reliably approximate the width of the 
zero mode band by: 

T = 2V2ke~ 2x (69) 
while the separation of the zeroth and the first band is 

ATi = ^ke~ x (70) 

with £1 the 1st root of J\{x). This characteristic behaviour of the above widths with respect 
to x implies that we can have a viable multigravity scenario if Tq 1 >, 10 26 cm and at the same 
time AIY 1 <, 1mm. This can happen for k ~ Mpi and x ~ 68. 

The band structure has the following form. The first bands have very narrow widths Tj 
with spacing ATi between them (i.e. between the z-th and the [i — l)-th band), which are 
approximately: 

Ti » e^ke~ 3x (71) 
« fa-^ke-* (72) 

where £ 2 n+i is the (n+l)-th root of J\(x), £ 2n is the -th root of ^2(^)1 an d e 2n+ i = j j^llr+i) ' 
e 2n — 7 ^ 2 |! 2 "v As we move on to higher bands their width is increasing and the spacing 
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between them is decreasing. In the limit m —>■ oo, the spacing disappears, i.e. ATj — ► 0, 
while their width tends to r, — > ^ke~ x . In this limit the function /(m) in (|68|) tends to 
/(m) = cos(2;?/m) and we have no forbidden zones. 

In order to find the gravitational potential we should at first find how the KK modes cou- 
ple to matter and determine their spectral density. Instead of working with wavefunctions 
in an infinite extra dimension with continuous normalization, it is more convenient effec- 
tively to compactify the system to a circle of length 2A/7, by assuming that we have a finite 
crystal of Af cells [2jl,2(J + 1)1]. We shall find the gravitational potential is independent 
of Af, so we will be able to take the limit Af — > oo and decompactify the geometry. 

In this regularization scheme we can readily see that the ever-present zero mode, corre- 
sponding to the 4D graviton, has the wavefunction: 

*<°-*>=eS (73) 

and its coupling to matter on the central positive brane is: 



The wavefunction eq([75|) is a smooth limit of eq(|5^) for m — > 0. Indeed, for large x we 
find the normalization constant in (p4[) in the region of the zero mode band is approximately: 

= fnTW (75) 

The approximation of the Hankel functions with their first term of their power series proves 
our previous claim. Furthermore, we can see that the KK states in the whole zero mode 
band for large x will couple to matter in the central positive brane as: 

a(m) = 7WW/ t{m> < 76 > 



with <f>(m) an unimportant m-dependent phase that doesn't appear in observable quantities. 

In the limit m — » oo the Hankel functions can be approximated by their asymptotic 
form. The wavefunction is then almost constant in z-space and the KK states couple to 
matter on the central positive brane as: 

aim) = -^^=L — e^ (m) (77) 



where again <fi(m) is an unimportant m-dependent phase. The couplings in the low lying 
bands interpolate between the above limiting cases. 
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Figure 6: Spectral density for the first three bands 



The standard procedure to calculate the spectral density is at first to find the number 
of states with masses smaller than m and then differentiate with respect to m. In our 
regularization prescription it is straightforward to count the states in the quasi-momentum 
space with masses lighter that m. The latter is simply: 



1% 



tot 



m) 



so the spectral density is its derivative: 

dN tot {m) Mzi dq 



p(m) 



dm 



71 



(=F) 



-q 



df(m) 



7T dm 2n y/\ - f{m) 



dm 



(7E 



(79) 



where the plus and minus signs will succeed each other for neighboring bands, starting with 
minus for the zero mode band. For the zero band we can have a reliable approximation for 
large x of the above function: 

1 1 



p(m) = J\f 



7lT 



(m/T f 



0) 



Thus the spectral density diverges as m — > T but the divergence is integrable and 
doesn't cause any problem to the following calculations. Actually, the spectral density 



diverges at the edges of every band since at these points ^/l — f(m) 2 
The only point that this doesn't happen is at m = where — > 
the result is finite. 



while 



df(m) 
dm 



^0. 



in such a way that 



In the limit m 

density: 



oo the bands disappear and the spectrum has constant spectral 



p(m) = Af 
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The gravitational potential, taking into account the Yukawa suppressions of the KK 
states, is simply: 

poo —mr 

V(r) = - dm\a{m)\ 2 p{m) (82) 

Jo r 

For distances r 3> AT^ 1 we have the effective potential: 

t 7-/ \ 11 f To , pim) _ mr 1 1 f 1 

M*rJ M ttM # 2 r J ^/l^ 

\l Q {Y Q r) - L Q (T Q r)\ (83) 



2Af 2 



r 



where Jo is the 0-th modified Bessel function and L is the -th modified Struve function. 
The gravitational potential is M -independent as expected. For distances AT^ 1 C r C T 1 
the above functions tend to I (T r) — > 1 and L (r r) — > 0, so we recover the 4D Newton 
law with G N = : 

V(r « To 1 ) = (84) 

On the other hand, for distances r 3> Tq 1 the asymptotic expansion of the difference of 
the modified Bessel and Struve functions gives a 5D Newton law: 

^» r »')-f^ < 85 > 

In case we take the crystal to be finite, there will appear a region for r ^> m\ x where 
gravity will turn again to 4D. However, as explained in the previous section, this region is 
well above the universe horizon and thus of no phenomenological interest. 

Finally, for distances r <, Ar^f 1 we will start to feel the short distance Yukawa type 
modifications to gravity due to the presence of the bands above the zeroth one. As in the 

zero hierarchy "bigravity" and " -\ h" "multigravity" model we expect these corrections 

to be important at scales of order the Planck length. 



5 Discussion and Conclusions 



We have shown that in models in which negative tension branes are sandwiched between 
positive tension branes there are anomalously light KK states that can lead to modifications 
of gravity. So far we have implicitly assumed that the gravitational force generated by a 
massive graviton is identical to that generated by a massless one but this is not, in general, 
the case because the massive graviton has additional degrees of freedom which do not 
decouple in the massless limit. As a result the interaction generated by a massive graviton 
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violates the normal 4D relation between the gravitational interactions of matter and light. 
This may be seen explicitly from the form of the massive and massless graviton propagators. 
Up to terms involving momentum vectors which do not contribute to T^G^^T^ due to 
momentum conservation, the propagator has the form [pi], |19|] 



1 ] J (2tt) 4 p 2 -m 2 -ie 



where for m ^ 0, t — 1/3 but for m = 0, t = 1/2. The difference between the two 
propagators 

SGT^ix - x') = - f dp * J"" 9 ? e -^-^) 



6 J (27r) 4 p 2 — m 2 — ie 
is due to the additional helicity components needed for a massive graviton. There are three 
additional components needed. Two, corresponding to the graviphotons, decouple at low 
energies as they are derivatively coupled to the conserved energy-momentum tensor. The 
third, corresponding to a scalar component, does not decouple and is responsible for 5G. The 
contribution of 5G to the one-graviton exchange amplitude for any two four-dimensional 
sources and T 2 ^ vanishes for light for which T = = but not for matter because 
T = Tq ^ 0, hence the normal 4D gravity relation between the interactions of matter 



and light is violated as was first pointed out by Dvali et al |[20|| . This is inconsistent with 
observation, for example the bending of light by the sun. 

It has been argued that this violation does not occur in the models discussed above. Two 



mechanisms have been identified. The first [21] involves a cancellation of the contribution 



from SG due to an additional scalar field, the radion, the moduli field associated with the 
size of the new dimension. This cancellation was related to the bending of the brane due to 
the matter sources following the approach suggested by Garriga and Tanaka p2" |. This field 



must have unusual properties because it is known that a normal scalar contributes with the 
same sign as SG and therefore cannot cancel it. An explicit computation of Pilo et al [^3| 
of the properties of the radion in the GRS model show that it indeed has a negative kinetic 
energy term leading to a contribution which cancels the troublesome 5G contribution. As 
a result the theory indeed gives rise to 4D gravity at intermediate distances. The presence 
of the ghost radion field is worrying for the consistency of the theory for it gives rise to 
"antigravity" at extremely large scales when the theory become five dimensional. Moreover 
as the ghost energy is unbounded from below it is likely that any theory is unstable when 
coupled to such a state [Q. It has been suggested that the problems associated with the 



ghost apply to all theories in which there is a negative tension brane sandwiched between 



positive tension branes because such theories violate the weak energy condition []i7 , [25 



The latter has been shown [^6| to require o~"(y) > and, as we have seen, models which 
have intermediate negative tension branes violate this condition. However it is not clear 
that theories that violate the weak energy condition are unacceptable p8[ . Indeed it has 
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been demonstrated that any theory involving a scalar field necessarily violates the condition 
p9| and the condition also does not apply to the gravitational sector |30|. Certainly more 



work is needed to clarify these issues. 



The second mechanism capable of cancelling the SG contribution was identified in [|TT 
offers some hope that the problems associated with the ghost state may be avoided. The 
mechanism follows from the underlying 5D structure of the theory. The graviton propagator 



in five dimensions has the tensor structure 19 



qMN,PQ „ [ 1/2 ( g MP g NQ + g MQ g NPj _ fg MN g PQ^ + (pp/p 2 ) (87) 

where t = 1/3 i.e. the same as for the massive four-dimensional graviton. The full five- 
dimensional one graviton exchange amplitude has the form 



T l MN G MN ^Tl Q oc Tl 1N [1/2 (g MP g^ + g^g NP ) - tg MN g^} T? 



PQ 



where Tmn is the five dimensional stress-energy tensor. In the compactifications considered 
here T^ 5 vanishes but in general T 55 does not. Thus the total amplitude can be written as 



T^T 2 ^ - tT l T 2 + (1 - t)T 5 15 T 5 25 - t(T 5 15 T 2 + T l T 2b ) 



ftu - 2 



+ (1/2 - t)T x T 2 + (1 - t)T 5 15 T 25 - t(T 5 15 T 2 + T X T 5 25 ) (89) 

The first term in square brackets is the normal four- dimensional amplitude corresponding 
to massless graviton exchange. The second term takes account of the 5G contribution as 
well as the contribution from T|. It may be written in the form 

i(T 1 -2T 5 15 )(T 2 -2T 25 ) (90) 



In [n| it was pointed out that this contribution necessarily vanishes for stable brane config- 
urations. To see this we first note that all observable amplitudes involving states confined 
to the visible sector brane have zero momentum transfer in the fifth direction. This means 



that we have to integrate over the fifth coordinate / dy V-G (5) (T - 2T|). However it 



was shown in [g7| that this combination vanishes in any theory which stabilises the extra 
dimension. As a result the contribution of the second line of eq(^) corresponding to SG 
vanishes and one is left with the usual four dimensional gravitational interaction. Note 
that this explanation does not require a field with negative kinetic energy as it arranges for 
the additional (scalar) graviton component associated with a massive graviton to decouple 
from matter via a cancellation of the various contributions to the stress energy tensor. 

As discussed above there may still be a troublesome negative kinetic energy radion 
field but now it too decouples because it also couples to f dy V—G^ (T — 2T| ). This may 
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eliminate the problems associated with a negative kinetic energy scalar. Indeed it suggests 
that we should be able to introduce an additional field to cancel the ghost field without 
disturbing the 4D structure of the gravitational interactions. 

In summary, in this paper we discussed three models based on different brane configura- 
tions, the three-brane " + H — " model, the four-brane " H \-" "multigravity" model and 

the crystalline "multigravity" model. The " + H — "model provides an instructive example 
showing that we cannot have special light KK states (and thus "multigravity") without 
intermediate "— " branes. This model, however, does provide a way of having the visible 
sector on a "+" brane with an hierarchical "warp" factor. The second " -\ \-" "multigrav- 
ity" model provides an interpolation between the KMPRS and the GRS models and leads 
to a new possibility in which gravity changes from 4D to 5D and then back again to 4D. 
The third model, the infinite crystalline "multigravity" model provides a further example 
of how the "multigravity" scenario appears every time that we have "— " branes between 
"+" branes. Again gravity changes from 4D to 5D at ultralarge distances (and back to 4D 
if the crystal is finite). While models with intermediate negative tension branes have very 
interesting phenomenological implications there is a question whether they are consistent 



as they necessarily violate the weak energy condition p5], |26fl . Related to this is the fact 
that they may have ghost states with negative kinetic energy. We have shown that if the 
brane configurations are stabilised these ghost states decouple, offering the possibility that 
the problems associated with such states may be avoided. Moreover the same decoupling 
condition ensures that the gravitational interactions associated with the KK excitations of 
gravity responsible for "multigravity" couple in the usual 4D way to matter and radiation. 
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